Large Non-Gaussianities with Intermediate Shapes from Quasi-Single Field Inflation 
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We study the slow-roll inflation models, where the inflaton slow-rolls along a trajectory whose 
orthogonal directions are lifted by potentials with masses of order the Hubble parameter. In these 
models large non-Gaussianities can be generated through the transformation from the isocurvature 
modes to the curvature mode, once the inflaton trajectory turns. We find large bispectra with a one- 
parameter family of novel shapes, interpolating between the equilateral and local shape. According 
to the in-in formalism, the shapes of these non-Gaussianities are different from a simple projection 
from the isocurvature non-Gaussian correlation functions. 



INTRODUCTION 

Inflation jl) is the leading candidate for creating the 
homogeneous and isotropic universe and generating the 
primordial density fluctuations for large scale structure. 
The condition for inflation to happen is that the inflatons 
stay near the top of the potential for a sufficiently long 
time, so that the vacuum energy drives the accelerated 
expansion of the Universe. 

Generically, it is found that such a condition needs to 
be fined-tuned, at least for various types of models that 
have sensible UV completion [|| . For slow-roll inflation, 
this means that the generic inflaton potential has a steep 
shape, characterized by a mass of order the Hubble pa- 
rameter H . On the other hand, it is general that the 
inflaton can roll in a space with multiple light fields. So 
at least one of the directions needs to be fine-tuned. A 
combination of these two aspects suggests a generic sit- 
uation - the inflaton slow-rolls along a trajectory whose 
orthogonal directions are lifted by potentials with masses 
of order H . We call this class of models the quasi-single 
field inflation. If the inflaton trajectory is straight, this 
is equivalent to the single field inflation, which gener- 
ates unobservably small primordial non-Gaussianities [[|. 
However, as we will show in this paper, once the trajec- 
tory turns, large non-Gaussianities with novel shapes can 
be generated, which are potentially detectable in current 
and future experiments. 

We call the mode in the tangential direction the curva- 
ture mode, and the modes in the perpendicular directions 
the isocurvature modes. The non-Gaussianities origi- 
nated from the curvature mode are very small because 
the potential is flat. Since the isocurvature modes in our 
model are massive, m ~ H, and can have large higher or- 
der interactions, the non-Gaussian correlation functions 
of the isocurvature modes can be very large. However, in 
dS space the amplitude of the quantum fluctuations of a 
massive field decays exponentially after horizon exit, and 
is also fast oscillating if m 3> H . This is why such modes 
are usually not considered. But for m ~ H , the sup- 
pression due to oscillation is only marginal. In the mean 



while, although the isocurvature amplitude still decays, 
the part of it that is transferred to the tangential direc- 
tion through turning trajectory becomes part of the cur- 
vature mode, which remains constant after horizon exit. 
The non-Gaussianities in the isocurvature modes are thus 
transferred in this way. This is the main physical picture 
behind our calculation. 

We use the in-in formalism 0] to compute the precise 
momentum dependence (shapes) of the bispectra. In this 
formalism, the natural way to implement the transforma- 
tion from the isocurvature to curvature mode is to intro- 
duce a two-point transfer vertex between the two modes. 
It is then straightforward to do a perturbative calculation 
for correlation functions according to the usual Feynman 
diagrams. As we will see, the effect of the transfer vertex 
is not a simple projection of the three-point function of 
the isocurvature modes. 

During the computation, we find that two equivalent 
representations of the in-in formalism, namely the origi- 
nal factorized form and the commutator form, are com- 
putationally advantageous in complementary ways. In 
certain parameter space, each representation encounters 
spurious divergence either in IR or UV. These divergences 
are called spurious because they will eventually be can- 
celled, but significantly complicate the analytical and nu- 
merical calculations. They are completely absent, how- 
ever, when viewed in a different representation. 

We find a one-parameter family of bispectrum shapes, 
lying in-between the well-known equilateral and local 
shape. The shape sensitively depends on the mass of the 
isocurvature mode; the size depends on the turning an- 
gular velocity and the cubic interaction strength among 
isocurvature modes. 

The mechanism that the isocurvature modes, in par- 
ticular their non-Gaussianities, can be transferred to 
the curvature mode through turning trajectory has been 
pointed out and studied intensively in the past |^-|^] • The 
focus was on the light isocurvature modes with mass 
much less than H . These modes do not decay after 
horizon-exit and non-Gaussianities come from non-linear 
classical evolution of super-horizon modes in multifield 



space. The resulting shapes of non-Gaussianities are 
local. Also note that for multifield slow-roll inflation 
models, it is found to be very difficult to generate large 
non-Gaussianities through this mechanism because the 
restrictive slow-roll conditions in all directions. Here in 
the quasi-single field inflation models, we focus on mas- 
sive isocurvature fields, which, as we emphasized and will 
show in details later, are crucial for generating large non- 
Gaussianities with new intermediate shapes. 



THE MODEL, MODE FUNCTIONS AND 
TRANSFER VERTEX 



-X- 
(a) 



5<T 

-x — X- 



(b) 



So i 
So , X ^Sa 



(e) 



FIG. 1: Feynman diagrams for the transfer vertex (a), and the 
corrections from the isocurvature mode to the power spectrum 
(b) and bispectrum (c). 



We consider a two-field model as an example. It is con- 
venient to write the action in terms of the fields in polar 
coordinates, 9 and a, that are tangential and orthogonal 
to the turning trajectory, respectively, 



S 



d^Xs 



-V sr {6) - V(a)} , (1) 



where R is the radius of the turning trajectory. The 
V sr (9) is a usual slow-roll potential and we choose the 
rolling velocity 9q > 0. The potential 



V(a) = \mW 



-V"'c 
6 



(2) 



traps the a field at a = |T(J. In principle, the param- 
eters that characterize the turning trajectory, 9q and R, 
and the cubic interaction strength V'", can vary along 
the trajectory. In this paper we consider the constant 
turn case in which they are all constant. 

We perform the perturbation in the gauge where the 
scale factor a(t) is homogeneous. The leading kinematic 
Hamiltonian density for the quantum fluctuations, 59 / 
and 8a i , in the interaction picture is 



Hn — 
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The leading interaction Hamiltonian density is 

Hi = -c 2 a 3 8a I 89 I , 
Hi = c 3 a 3 8a 3 j , 

where c 2 = 2R9 , c 3 = V"'/6 are constants. 

We quantize the fields 89i and 8oi in the momentum 



.(3) 



(4) 
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space p, 



+ u*_. 
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K = 



6 f 
-p -p' 



where a p and b p are independent of each other and each 
satisfies the usual commutation relation. The mode func- 
tions in terms of the conformal time r = J dt/a(t) are 



-)^H^(-pr) , (7) 




where v = m 2 / H 2 and m 2 = m 2 + 7 9%. As men- 

tioned in the Introduction, the amplitude for the isocur- 
vature mode v p decays as (— r) 3 / 2_!y at late time r — » 0. 
If m 2 /H 2 > 9/4, v is imaginary, v p has an additional fast 
oscillating factor, ~ e I/ln (-P T / 2 ) ; even after the horizon 
exit. This suppresses its contribution. We will consider 
the case < v < 3/2, i.e. 9/4 > m 2 /H 2 > 0, in this 
paper. 

The term that is responsible for the transformation 
from the isocurvature to curvature mode is (||). This in- 
troduces the "transfer vertex" (Fig. 0(a)). The contribu- 
tion from isocurvature to curvature correlation functions 
can therefore be computed according to the Feynman 
diagrams in Fig. (b)(c). The term (§) is the leading 
source for large non-Gaussianities, since all other cubic 
interactions in the expansion involve 89 which is in the 
slow-roll direction. 



BISPECTRA 

We compute the three-point correlation function (89 3 ). 
The bispectrum (( 3 ) is related to it through the usual 
time-delay relation £ w —H89/9. 

The in-in formalism [H gives 



89 3 {t) 
10) , (8) 



Texp ( i / dt'H^t') 
7>xp j / / dt'Hi(t' 

to 



where t is the end of the inflation, and t is the infinite 
past. The Hj is the interaction Hamiltonian, 



H I= d 3 x (Hi+Hl) 



(9) 



One can directly expand the exponentials in (@). For our 
purpose, the relevant terms are those that involve four 
ii/'s. Each resulting quartic integral contains two "fac- 
tors" . Each factor has an (anti-)time-ordered integration 
but there is no cross time-ordering between the two. We 
call this form the "factorized form" . One can also rewrite 
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all the integrands into one single time-ordered quartic in- 
tegral, 

dt 2 I dtz I dt^ 

([Hi (U)][Hj (t 8 ),Vi (i 2 )° [Hj (h ) , 60! (t) 3 ]]]]) . (10) 

We call this the "commutator form" . The Feynman dia- 
gram Fig. 0(c) corresponds to replacing one of four i/j's 
in (^) or QIC] ) with H% and the rest three with H\. A 
straightforward exercise of Wick contraction leads to, for 
the factorized form, 

- l2clc 3 u* pi u P2 u P3 (Q) 



x Re 



dfi a v pi u pi (fi) / df 2 a v Pl v P2 v P3 (f2) 



I dr x a 3 v* 2 u'* 2 (n) / dr 2 a 3 v* 3 u'* 3 (r 2 ) 

) J —oo 

x (27r) 3 <5 3 (^ pi) + 9 other similar terms 
+ 5 permutations of Pi ; (11) 
and for the commutator form, 
12clc 3 u Pl (0)u P2 {0)u P3 (0) 

/0 /•Tl />T 2 /-T3 4 

dn / dr 2 / dr 3 / dr 4 Y[(a 3 {n)) 
-oo — oo J — OO J — OO • i 

x o(r 2 ) (u pi (n) ~ cc.) (^(1-1)1^(7-2) - c.c.) 

x (v P3 (t 2 )u* 3 (r 4 )w p * (t 4 ) - c.c.) t; P2 (t 2 )u* 2 (r 3 )ti£ (r 3 )] 

x (27r) 3 (S 3 (^ p,^ + 2 other similar terms 



5 permutations of p; 



(12) 



We write the argument Ti only once if they are all the 
same in one integrand, and the prime denotes the deriva- 
tive respective to the conformal time Tj. The full details 
are presented in Ref . [jy] . 

It is subtle to evaluate these integrals. Let us first 
look at the factorized form. In the UV limit, 75 — > —00, 
the integrands are fast oscillating. For the Bunch-Davies 
vacuum, the convergence of the integration is achieved by 
slightly tilting the contour clockwise or counter-clockwise 
into the imaginary plane, Tj — > — oo(l±ie). In fact, if the 
integration is convergent at IR, Ti — > 0, it is numerically 
much more convenient to do a Wick rotation Tj — > iZi for 
the left factor and opposite for the right. The integration 
range for z% is from —00 to 0. Now it is explicit that the 
integrand of each factor decays exponentially at UV. This 
is the case for < v < 1/2. 

However, for 1/2 < v < 3/2, using the asymptotic be- 
havior of the mode functions, it is easy to see that each 
term in (|ll]) is IR divergent. Physically, larger v corre- 
sponds to smaller m for the isocurvature mode. So the 
mode decays slower. The conversion from isocurvature 
to curvature mode lasts longer after the horizon exit. As 



we will see, this causes slower convergence in the IR, but 
not divergence. 

Let us look at this in the commutator form (12). Due 



to the subtraction of the complex conjugate in various 
terms, the integrand decreases faster in IR. For example, 
u' Pi (ti) behaves as (— n) as t\ — > 0; but after subtracting 
off its complex conjugate, we have (— ti) 2 . The next two 
terms are slightly more complicated but similar. Since all 
these three terms are imaginary, v P2 v* 2 u'* 2 in the 4th line 
must be imaginary to make the whole integration real. 
For T2,3 — > 0, the leading term of v P2 v* 2 u'* 2 is real, so we 
need the subleading term in this limit. This also increases 
the IR convergence. Overall, it is not difficult to see 
that the IR convergence is achieved for all < v < 3/2. 
However, the Wick rotation no longer works in this form. 
The original integrands from the left and right factor 
have been multiplied together, so after Wick rotation, 
the exponential decay of some factors are cancelled by 
the exponential growth of the others. 

To summarize, for 1/2 < v < 3/2, the factorized form 
is well behaved in UV but encounters spurious divergence 
in IR in the intermediate step, while the commutator 
form is well behaved in IR but the UV convergence be- 
comes tricky Q. Therefore we have proved that the 
whole integral has no divergence. In fact, to see both 
types of convergence at once, we can choose a cutoff t c 
in the middle, and make the whole expression take the 
factorized form in the UV, Ti < r c , and commutator form 
in the IR, Ti > t c Jllj]. This provides an efficient way to 
numerically calculate the shapes of the bispectra. 



SQUEEZED LIMIT AND SHAPE ANSATZ 

We now look at the squeezed limit, p% <C pi = p 2 - 
In this limit, simple analytical expressions for shapes are 
possible. This can also help us construct simple analyt- 
ical ansatz for the full shape, to facilitate future data 
analyses. In this limit, both Eq. (pd|) and (03) become 



c|c 3 



1 



HR 6 |+f 

p 2 p2 p2 



s(v) , 



(13) 



where s[y) is a pi-independent numerical number, but 
involve complicated integrals jll| . The result is presented 
in Fig. I 

Some comments are in order, to explain the behavior of 
s(v) near v = and v = 3/2. We have approximated the 
asymptotic behavior of Hi, > (—paTi) in the small p 3 limit 
as —i{2 v T(v)/-K)(—pzTi)~ u . For very small v ~ 0, this 
requires pzjpx <§C e -1 ^. So pi needs to be increasingly 
small as v — > 0. Otherwise, if we fix a small ps, near 
v = 0wc should instead use ?(2/tt) ln(— P3n) as a better 
approximation. Therefore the rising behavior in Fig. |2| 
near v = does not mean that the non-Gaussianities are 



3 




FIG. 2: The numerical coefficient s(y) in the squeezed limit. 



blowing-up, rather signals the change of the shape to 



7/2 3/2 
Pi P2P 3 



(14) 



As v — > 3/2, m — > 0, the curvaton fluctuations do not 
decay, so its conversion to the curvature mode diverges 
in the constant turn case and an e-fold cutoff is neces- 
sary. Interestingly, in this limit, our shape approaches 
the local form, as expected from the multi-field models 
studied in Ref. But here the non-Gaussianities can 

be made very large by having a large V'" , since we are not 
restricted to the slow-roll conditions even in the massless 
isocurvaton limit. 

Combining ([[3]) and (Q), a good ansatz for the full 
shape can be taken as, up to an overall amplitude, 



(P1P2P3) 3/2 

(Pi +P2 +PzY /2± " V V(Pl +P2 +P3) 



nP2Pz 



(15) 



where N v is the Neumann Function. This ansatz gives 
an overall good match with the numerical results |ll| . 

Note that in the squeezed limit, this one-parameter 
family of shapes goes as p^ 3 ^ 2 ~ u . This interpolates be- 
tween the equilateral form, pg , and the local form, pg 3 
fL3f , so we call it the "intermediate form" . Two examples 
of the shape ansatz are shown in Fig. 0. 




horizon exit, we find that its amplitude is decaying and 
its shape goes as pg 2v in the squeezed limit. So at least 
in this model, the shape of the correlation function is 
changed during the transfer. It is important to study 
this aspect in other multifield models, such as u%- 



SIZE OF NON-GAUSSIANITIES 

The size /nl of a bispectrum is defined by taking the 
equilateral limit JL5|, 



(C 3 



^/^p c W* 3 (Ep*) 



(16) 



where P^ is the power spectrum. Using the relation £ = 
-H56/6 and P c = H* /(4ir 2 R 2 d 2 ) sa 6.1 x 10" 9 , we get 




(17) 



We investigate the order of magnitude of each factor. The 
a(v) should be evaluated numerically, similar to s(u), but 
now in the equilateral limit. For example it is 0(1) and 
positive near v = 0. If we require that, in V(o), the 
quadratic term dominates over the cubic interaction for 
a < H, so that we can trust the mode function, we need 
\V'"\/H < (m a /H) 2 ~ 0(1). The perturbative method 
we used gives restriction on the size of 9 /H because this 
parameter determines the strength of the transfer vertex. 
For example, the correction to the power spectrum can be 
simply calculated using Fig. [I] (b), 8Pq ~ (Qq/H) 2 Pq and 
it is scale-invariant, so for it to be perturbative we need 
(9q/H) 2 < 1. It is possible that the non-Gaussianity 
is larger if 8 /H ~ 1, but to trust the perturbative re- 
sults in this paper, 80/H < 0(1). Overall we see that 
\ffi L \ < 0(1O 4 ), and its sign is the opposite of V". It 
will be very interesting to constrain it using the current 
and future data |l(|. It is also interesting to study what 
the natural values for V'" and 9q are from a more funda- 
mental theory. 
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FIG. 3: Shapes of bispectra with intermediate form: 1) quasi- 
equilateral (y — 0.2), 2) quasi-local [y = 1). The amplitudes 
are normalized by a factor of (pip-zP'i) 2 to be dimensionless. 

For comparison, we look at the three-point function of 
the isocurvature modes, (S<r 3 ). Evaluating it after the 
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